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1 Introduction 

The past two decades have seen significant progress in the understanding of the correlation 
functions of local operators in spin-^ chains. This report is about the extention of recent results 
for the ground state correlators of the XXZ chain, surveyed below, to finite temperatures. 

The development was initiated with the derivation of a multiple integral formula for the 
density matrix of the XXZ chain by the Kyoto school [24-26] which relies on the bosonization 
of q-vertex operators and on the (/-Knizhnik-Zamolodchikov equation [18, 35]. An alternative 
derivation of the multiple integral formula was found in [29]. It is based on the algebraic Bethe 
ansatz and made it possible to include a longitudinal magnetic field. 

The multiple integral formulae, however, turned out to be numerically inefficient. They were 
hence not much used before it was realized [9] that they may be calculated by hand, at least in 
principle. This result generalized after many years Takahashi's curious formula [39] for the 
next-to-nearest neighbour correlator and inspired a series of works devoted to the explicit cal- 
culation of short-distance correlators in the XXX [10, 1 1, 15,32-34] and XXZ chains [27,28,40]. 
It further triggered a deep investigation into the mathematical structure of the inhomogeneous 
density matrix of the XXZ chain, which was started in [12-14] and still continues [2-7]. 

In [2] a minimal set of equations that determines the inhomogeneous density matrix was 
derived and was termed the reduced ^-Knizhnik-Zamolodchikov (rqKZ) equation. The rqKZ 
equation made it possible to prove that the correlation functions of the inhomogeneous XXX 
model depend on a single transcendental function which is basically the two-spinon scattering 
phase. This was generalized to the XXZ and XYZ models in [3,7], where further transcendental 
functions were needed. 

A new 'exponential form' of the density matrix was derived in [5] and [4] for which the 
homogeneous (physical) limit can be taken directly. The most recent papers [6, 8] aimed at 
understanding how the exponential formula works in the 'free fermion' XX limit. This led 
to a novel formulation also for generic q. A crucial tool was a disorder field acting on half 
of the infinite lattice with 'strength' a. It regularized the problem further and simplified the 
exponential formula in a way that the exponent depends only on a single transcendental function 
to and on special operators b and c resembling annihilation operators of (Dirac) fermions. 

From the above studies we observe the following. In the inhomogeneous case the multi- 
ple integrals reduce to polynomials in a small number of different single integrals related to 
the correlation functions of only nearest- neighbouring lattice sites. These constitute a set of 
transcendental functions which determine what we call the 'physical part' of the problem. The 
coefficients of the polynomials are rational functions of the inhomogeneity parameters. They 
are constructed from various L-operators related to the symmetry of the models and constitute 
the 'algebraic part'. We call such type of separation of the problem into a finite physical part and 
into an algebraic part 'factorization', since it can be traced back to the factorization of multiple 
integrals into single integrals. We believe that factorization is a general feature of integrable 
models (for a similar phenomenon in the form factors for the Ising model see [16]). 

A generalization of the integral formula for the density matrix of the XXZ chain to finite 
temperature and magnetic field was derived in [19, 21, 22] by combining the techniques devel- 
oped in [29] with the finite temperature formalism of [30, 31, 36-38]. Remarkably, the form of 
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the multiple integrals for the density matrix elements is the same in all known cases. The phys- 
ical parameters (temperature T , magnetic field h, chain length L) enter only indirectly through 
an auxiliary function which is defined as a solution of a non-linear integral equation. 

The auxiliary function enters into the multiple integrals as a weight function. This implies 
that the factorization technique developed for the ground state correlators in [9] does not work 
any longer. In our previous work [1] we nevertheless obtained a factorization of the correlation 
functions of up to three neighbouring sites in the XXX model at arbitrary T, h by implicit use 
of a certain integral equation. Comparing the factorized forms with the known results for the 
ground state we could conjecture an exponential formula for the special case of T > but h = 0. 
Surprisingly, the formula shares the same algebraic part with its T = counterpart; one only 
has to replace the transcendental function by its finite temperature generalization. The results 
easily translated into similar results for the ground state of the system of finite length [17]. 

In this work we extend our analysis to the periodic XXZ chain 



in the antiferromagnetic regime (7 > and A = ch(ri) > —1) and in the thermodynamic limit 
(L = 2N oo). We identify an appropriate set of basic functions describing the neighbour cor- 
relators in the inhomogeneous case. The algebraic part of the problem without magnetic field 
is neatly formulated in terms of the operators b and c as in the ground state case. The meaning 
of the disorder parameter a, necessary for the construction of these operators, is yet to be un- 
derstood for finite temperatures. It, however, naturally modifies one of our auxiliary functions, 
the density function G and allows us to reduce the number of basic functions characterizing the 
physical part from two to one. 

Still, we go one important step further. We extend our conjectured exponential formula for 
the (finite temperature) density matrix such as to include the magnetic field. At first sight, this 
may seem to require only trivial modifications, as the Hamiltonian commutes with the Zeeman 
term. The magnetic field, however, breaks the Uq{sl2) symmetry and, as far as the factorization 
of the integrals is concerned, brings about serious difficulties even for the ground state corre- 
lator problem. For this reason an essential modification of the operator in the exponent of our 
exponential formula is required which leads to novel formulae even in the zero temperature 
limit. The prescription is, however, remarkably simple. We have to add a term whose algebraic 
part is determined by a new operator H, such that the operator in the exponent is now a sum of 
two ingredients. One is formally identical to the operator already present at vanishing magnetic 
field, the other one is constructed from H (note that even the former part is not independent of 
the field; it includes transcendental functions which are even functions of h). 

We finally point out a simplification compared to the ground state case, particularly relevant 
at finite magnetic field. Although we are dealing with highly nontrivial functions, all correla- 
tion functions should simplify in the vicinity of T = oo. Thus, the high temperature expansion 
technique can be applied to the multiple integral formulae at T > as was shown in [41,42]. 
We use this in order to test our conjecture for the exponential form of the density matrix. 

Our paper is organized as follows. In section [21 we recall the definition of the density matrix 
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and the multiple integral formulae. In section [3] we describe the basic functions that determine 
the physical part of the correlation functions. Our main result is presented in section |4] (see 
eqs. (l36l)- (|38l) ). It is a conjectured exponential formula for the density matrix of the XXZ 
chain at finite temperature and magnetic field. Section [5] is devoted to the simplest examples 
of correlation functions, the cases of n = 1,2,3, for which we show novel explicit formulae. 
In section [6] we summarize and discuss our results. Appendix lAl contains the proofs of two 
formulae needed in the main body of the paper, appendix |B] a derivation of the factorized form 
of the density matrix for n = 2 directly from the double integrals, and appendix O a short 
description of the high-temperature expansion technique. 



2 Multiple integral representation of the density matrix 

Let us recall the definition of the density matrix of a chain segment of length n. We would like 
to take into account a longitudinal magnetic field h which couples to the conserved z-component 

1 ^ 

of the total spin. Then the statistical operator of the equilibrium system at temperature T is 
given by 

e T ^ T 

PN{m = (3) 

tr-A^+i,...,A/ e T + T 

From this operator we obtain the density matrix of a chain segment of length n by tracing out 
the complementary degrees of freedom, 

D„(r,/z|A^) = tr_A,+i,..._Lo,«+i,n+2,...,A'Piv(7^,/2), n=l,...,N. (4) 

The density matrix Dn{T,h \N) encodes the complete equilibrium information about the segment 
consisting of sites 1 , . . . , n which means that every operator acting non-trivially at most on 
sites 1 , . . . , n has thermal expectation value 

{0)T,h = trL...,n{DniT,h\N)0). (5) 

We know a multiple integral representation for the density matrix @i in two limiting cases, 
the thermodynamic limit ^ oo [19, 21] and the zero temperature and zero magnetic field 
limit [17]. For the two limits we shall employ the notation 

Dn{T,h) = lim D„{T,h\N) , Dn{N) = lim \imDniT,h\N) . (6) 

These two density matrices are conveniently described in terms of the canonical basis of endo- 
morphisms on (C^)®" locally given by 2 x 2 matrices e^, a, P = ±, with a single non-zero entry 
at the intersection of row P and column a, 

D„{T,h) = Dn^\';:±{T,h) eil\ . ..e„il , Dn{N) = Dn^-±{N) eii\ . ..ent , (7) 
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where we assume implicit summation over all (Xj.^k = ±- We further regularize the density 
matrices by introducing a set of parameters ?ii , . . . , a in such a way that 



A. (r, h) = lim lim D„ {lu...X\T.h;a)eil\...ent, 

/.i,...,A„^0 a^O 

Dn{N)= lim lim3Pi;;;i(Xi,...,X„|A^;a)eiP^^^ 



(8a) 
(8b) 



From here on we shall concentrate on the temperature case (fSal ). Later we will indicate the 

modifications necessary for (|8bl) . We call D„ aj' ....a", (^i , • • • An I oc) the inhomogeneous den- 
sity matrix element with inhomogeneity parameters Xj. For a = it has a clear interpretation 
in terms of the six- vertex model with spectral parameters . . . on n consecutive vertical 
lines [22]. For h,T = the variable a can be interpreted as a disorder parameter [25]. In 
the general case we simply define the inhomogeneous density matrix element by the following 
multiple integral, 

s 

n 



Xj-1 



dco e""^ 



Xj-l 



det[-G(a)j,?iyt;oc)] 



ni<;<fc<,, sh(Xfc - ?i^) sh(a)^- - (Ofe - ri) 



(9) 



Here s is the number of plus signs in the sequence (aj)"^j, and s' is the number of minus signs 
in the sequence ((3j)"^j. The factor ds^m^s' reflects the conservation of the z-component of the 
total spin. For j = 1, ... ,5 the variable xj denotes the position of the jth plus sign in (aj)"^j 
counted from the right. For j = s+l,...,nit denotes the position of ( j — 5)th minus sign in 
((3j)"^j. The integration contour depends on r\. We show it in figureHJ This contour will also 
appear in the integral equations which determine the transcendental functions a, a and G and in 
the definition of the special functions in the next section that determine the physical part in the 
factorized form of the correlation functions. For this reason we call it the canonical contour. 
The integral equation for a is non-linear, 



Ina(X) = 



27sh2(ri) 



d(0 sh(2ri)ln(l + a(oo)) 



rsh(X)sh(?i + ri) Je Ini sh(?i-a) + ri) sh(X-(o-ri) 



(10) 



There is a similar integral equation for a (see [20]), however, since a = 1 /a we do not need to 
consider it here, a is usually called the auxiliary function. The combination 1 / ( 1 + o) has a 
natural interpretation as a generalization of the fermi function to the interacting case [23]. Note 
that the right hand side of equation (flOl) is the only place where the thermodynamic variables T 
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Figure 1 : The canonical contour C for the off-critical regime A > 1 (left) and for the critical 
regime — 1 < A < 1 (right). 



and h enter explicitly into our formulae for correlation functions. They neither enter explicitly 
into the multiple integral formula ^ nor into the linear integral equation for G which is 

G(X,/y;a) = -coth(X-A*)+e«^coth(^-A/-ri)+ / p^^^il^Kil-(o;a) . (11) 

G can be interpreted as a generalized magnetization density (see [20]). Compared to our previ- 
ous definition [20] we introduced the additional parameter a here which also enters the kernel, 

K{X- a) = e«^ coth(X - ri) - e""^ coth(X + ri) . (12) 

An equivalent integral equation for G which uses a instead of a and which is sometimes useful 
is 

G(X,^;a) = -coth(X-^)+e-«^coth(X-^ + ri)- / ^^^^:^K{X-^-a). (13) 

Setting a = the function G(k,ij; a) turns into the function G(X,;u) which played a crucial role 
in our previous studies [1, 20, 22]. We have introduced a in such a way into ^ and ([TTI) that 
for T,h = the multiple integral representation Q turns into the finite-a expression that can 
be obtained within the ^-vertex operator approach of [25]. Our main motivation for introducing 
a into our functions was to enforce compatibility with the formalism developed in [6], where 
a is an important regularization parameter. The usefulness of this modification will become 
clear in section HI The parameter a will allow us to write our formula for the density matrix in 
factorized form in a very compact way. 

Let us briefly indicate the changes that are necessary in the finite length case (|8b| ). It turns 

out [17] that Aia'i' (^1' • ■ ■ A«|^;Oc) has a multiple integral representation of the same form 
as (|9l), that even the integral equation for G remains the same and that the only necessary mod- 
ification is in the driving term of the non-linear integral equation ([TOl) , where the physical pa- 
rameters enter, which in this case are the length L = 2N of the chain and an arbitrary twist 
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4> G [0, 27r) of the periodic boundary conditions (for details see [17]). The non-linear integral 
equation for the finite length case is 

, o-^ r /■ d« sh(2ri)ln(l + a(a))) 

lna(X) = -2i4> + Lr|+Lln — ^- — ^ - / — -rr^^ — ^ ' , . r. (14) 

^ ' ' Vsh(X+5)y ye2rash(?i-a) + ri)sh(X-a)-ri) 

When we derived the multiple integral representation ^ in [19] and [17] we assumed that 
the inhomogeneity parameters Xj are located inside the integration contour C. This has to be 
taken into account when calculating the homogeneous limit in (|8b1) . where the canonical contour 
should be first shifted to ±ri/2. 



3 The basic functions 

In this section we describe the functions constituting the 'physical part' of the factorized cor- 
relation functions of the XXZ chain at finite T and h. A description of the algebraic part will 
be given in the next section. According to our experience the physical part of the correlation 
functions can be characterized completely by two transcendental functions cp and (O. 
Let us start with the more simple function 

(pG";cx = 1+ / 7,, ■ (15) 

^^"^ ' ie 7rz(l + a((o)) 

This function is related to the magnetization m{T,h) through (p(0;0) = —2m(T,h) which we 
expect to belong to the physical part if the magnetic field is non-zero. 
In order to introduce the function (0 we first of all define 

V(^i,^2;a) = |^^^^^^(-coth(o)-^2)+e-«^coth((o-^2-il)). (16) 

Those readers who are familiar with our previous work [1] will recognize this as the anisotropic 
and 'a-deformed' version of the function \|/(/ji,/J2) introduced there. The function (0 is a modi- 
fication of \|/ obtained by adding and multiplying some explicit functions, 

co(^i,^2;a) = -e"^^i"^2)^^^j^^2;a) - - — -j-^K{^i -^2\-a) . (17) 

2 cosh (^) 

Here K{X;a) is the kernel defined in ([T2)) . The relation between a)(^i,/j2;a) and \|/(/ji,/j2;«) 
is similar to the relation between 7(^1,^2) and \|/(/ji,/J2) in the isotropic case [1]. The function 
CO is closely related to the neighbour correlators (see appendix |B]). In the critical regime for 
r, — > it becomes the function (o(^, a) of the paper [6] if we set ^ = e^'^^^^. 
An important property which follows from the definitions ([T2|) and ([T6l) is that 



03(^/2, A^i; -a) =a)(;Ui,;U2;a). (18) 



It implies 

a)(^2,A'i;0) = (o(a^i,a'2;0) , d{^i2,ni;0) = -(o'(a/i,a'2;0) , 



(19) 
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where for later convenience we introduced the somewhat unusual notation 

a)'(A/i,A'2;a) =aa(e'^(^2-^")a)(A'i,A'2;a)) . (20) 

At this point we would like to stress that the physical parameters T , h or A^, respectively, do 
not enter the definitions of cp and (O explicitly. The basic functions defined in this section are 
therefore suitable for both, the finite temperature and the finite length case, the only distinction 
being the use of different auxiliary function (fTOl ) and (fT4l) . respectively. 

In the high-temperature limit (see appendix O we observe that 

ga(^i-^2) T /an\ / 1 \ 

co(A/i,A/2;a) = tanh2(^yjiS:(^i-^2;-a) + Of -j . (21) 

Using eq. (|2T)) we conclude that both functions (o(;Ui,/J2;0) and (jo'(/ji,;U2;0) do not have zeroth 
order terms in their high-temperature expansions 

a)(^i,A/2;0) = o(i/r) , (o'(a/i,a^2;0) = o(i/r) . (22) 

The same is true for the function cp, 

(p(;u;a) = 0(l/r). (23) 

We mention the properties of these functions for a = with respect to reversal of the mag- 
netic field; cp(A';0) is an odd function of h, \|/(/ji,/J2;0) and 3aV(A'i7A'2;oc)|a=o are even. These 
properties will be implicitly used below. The proof relies on the simple fact that the quantum 
transfer matrix (or its slight generalization, see below) associated to the present model respects 
the spin reversal symmetry, and therefore the eigenvalues are even functions of h. 

Once this is realized, the proof for cp(/j;0) is rather obvious. One only has to remember the 
relation between cp(/j;0) and the largest eigenvalue A(/j) of the quantum transfer matrix, 

(p(A/;0) = r^lnA(^). (24) 

The above argument then implies that cp(/j;0) is odd with respect to h. 

The proof for \\f{iJi,iJ2;0) is less obvious. We first of all introduce a generalized system. 
Consider an 'alternating' inhomogeneous transfer matrix. In the framework of the quantum 
transfer matrix, we associate spectral parameters in alternating manner {u, —u, u,—u---) to 2X 
vertical bonds, while keeping the spectral parameter on the horizontal axis fixed as ^2- Next 
we add 2M vertical bonds and associate with them spectral parameters again in alternating 
manner, {u' + ^\,^\ — u' ,u' -\- ^\,^\ — u' ■). We then take the limit X, M — oo under the 
fine tuning, luJi = 2|37shri, 2a'M = — 257shri. Note that the original system is recovered by 
taking 5 = 0. By neglecting the term depending on the overall normalization, one obtains the 
following expression for the modified largest eigenvalue A(^2,j"i) of the generalized quantum 
transfer matrix. 
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The modified auxiliary functions a((0,/Ji), a((0,/Ji) satisfy equations similar to (fTOl) . and the 
equation for the latter is relevant here, 

lna(X,A/i) = ^-^^^^^e(^)+257sh(ri)e(^-A/i)+ / ^if(^-a);0)ln(l + a(a),;/i)) . 
11 Je 2ni 

(26) 

We take the derivative of both sides of (|26l ) with respect to 5, 

a(X,A/i)=27sh(Ti).(X-;.i)+ / ^^(^-«;0)7l4^^77^. (27) 

J e 2711 (1 + a(a),Aii)) 

where a(?i,^i) := 5(a.^^,) One compares (fTTI) with (|26|) and concludes 

a(;i,A/i)=27sh(ri)G(X,A/i;0). (28) 

Similarly we take the derivative of lnA(//2,/^i) with respect to 6 and find 

5, w X /■ d(o a((0,A/i) /■ dco , , 2ysh(ri)G((D,A/i;0) 

^lnA(u2,A/i) =- / — e((£)-u2)- / , = - / — efw-m) — , / 

= -7sh(Ti) / ^ (?^!!':uf\^ (coth(a)-,/2-Tl) -coth(a)-^2)) , (29) 

where we have used (l28l) in the second equality. By comparing the above equation with (fT6l) . 
one obtains 

V(a'i,A'2;0) = j-j^^lnA(A/2,A'i)l5=o- (30) 

Then the evenness of \|/(/Ji,/j2;0) follows from the same property of the generalized transfer 
matrix. 

Finally we show that , X2; cx) |a=o is also even. To prove this we consider the relation 

(IB. 51) in appendixlB The Ihs, D^^(?ii,?i2) +7)^i(^iA2), is invariant under + ^ — , hence it 
is even with respect to h. The first term in the rhs is also even as it is proportional to ,^2; 0) 
(see (IB.l II) ). Thus, the content of the bracket in the second term of the rhs should be also even. 
Thanks to (IB. 31) and (IB. 171) it is represented as 

D+(Xi)+D+(X2)-2D++(XiA2) = ^^^¥(^iA2;0) + '°'^^^^^ 
Thus, we conclude that 3aV|/(?iiA2;oc)|a=o is even. 



4 Thermal correlation functions of local operators 

In this section we are formulating our main result which is a conjectured explicit formula for 
the correlation functions of local operators in the XXZ chain at finite temperature and finite 
magnetic field. The sources of this conjecture are the results of the previous two sections that 
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followed from the finite temperature algebraic Bethe ansatz approach of [17, 20, 22] and the 
results of [4,6,7], where the exponential formula was discovered as a consequence of studying 
the rqKZ equation. Unfortunately, both approaches differ considerably in spirit and notation. 
We will try to reconcile them while keeping as much as possible of the original notation. We 
have to ask the reader to be forbearing though if this sometimes leads to confusion. 

In [6] much emphasis was laid on developing a formalism which applies directly to the 
infinite chain with lattice sites j E Z. To keep things closely parallel we therefore concentrate 
in this section on the temperature case and comment on the finite length case only later in 
section [6l All operators which act non-trivially on any finite number of lattice sites span a 
vector space W. Because of the translational invariance of the Hamiltonian we may content 
ourselves (as long as we keep a = 0) with operators which act non-trivially only on positive 
lattice sites, j E N. We shall denote the restriction of to the first n lattice sites by 0[i „]. The 
inhomogeneous density matrix satisfies the reduction identity 

tinDn{Xi,...,Xn\T,h;0) =D„_i(Xi,...,X„_i|r,/z;0). (31) 
It follows that the inductive limit 

lim tri,...,„ {D„{li ,...X\T,h;0)Oun]) (32) 
exists and defines an operator ^ : W ^ C such that 

D*T,hi^) = {0)T,h (33) 

is the thermal average at finite magnetic field of the local operator in the inhomogeneous 
XXZ model. Note that 

. ..e„l') = D,,l\'--±{Xu . . . X\T,h;0) . (34) 

For this reason we may interpret Dj as a kind of 'universal density matrix' of the XXZ chain. 
Let us define a linear functional tr : W ^ C by 

tr(0) = ...itri itr2 ^tr3...(0). (35) 

with try the usual traces of 2 x 2 matrices. Then we conjecture that an operator Q. -.W ^ W 
exists such that Dj ^ = tr e^. More precisely we propose the following 

Conjecture. For all eW the density matrix Dj ^ can be expressed as 

D^,,(0)=tr(e^(0)), (36) 

where tr is the trace functional d35l) and Q.:VQ is a linear operator that can be decom- 
posed as 

a = 0.1 + 0.2 (37) 
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with 

ai = -\im[ /^^a)(^i,^2;a)b(Ci;a-l)c(C2;a), (38a) 

^22 = -lim /^(p(A/i;a)H(^i;a). (38b) 

Here L,j = o^j, j = 1,2, and a)(^i,/J2;0^) (^^^d (p{pi;a) are the functions defined in 4771) and 
rfTTl) . The operators b, c and H Jo no? depend on T or h. They are purely algebraic. Their 
construction will be explained below. The integrals mean to take residues at the simple poles of 
b, c and H located at the inhomogeneities l^j (see below). 

In fact, the operators b and c are the same as in the ground state case [6]. The operator H is 
new in the present context but can be defined using the same algebraic notions underlying the 
construction of b and c. Note that lim/,^o ^{l^\ 0) = which implies that lim/,^o ^2 = 0. Hence, 
as in the isotropic case [1], we observe that the algebraic structure of the factorized form of the 
correlation functions is identical in the ground state and for finite temperature as long as the 
magnetic field vanishes. Due to the properties of the function (0 we recover the result of [6] 
in the zero temperature limit at vanishing magnetic field. In the high-temperature limit, on the 
other hand, we conclude with (l22l) . (1231) that lim^^oo f2 = and that all correlation functions 
trivialize in the expected way, 

limD*,, = tr. (39) 

For the definition of the operators b, c and H we first of all generalize the space of local 
operators W to a space of quasi-local operators of the form 

e«^lL-~^^0, (40) 

where is local, and denote this space by Wa. The operators b, c and H then act as 

b(C;a):Wa^W«+i, c(C;a):Wa^Wa_i, H(C;a):Wa^Wa (41) 

which implies in particular that b(^i;a — l)c(^2;oc) : W^. 

The z-component of the total spin is the formal series 5^ (see equation We denote its 
adjoint action by 

§(X) = [5i,X]. (42) 
Then q""^ : Wa. The spin reversal operator defined by 



7eZ 



(43) 



clearly is a map J : Wa — > W-a. 

The operators b, c and H will be defined in two steps. We first define endomorphisms 
h^ki], C[;t/] ^iid Hjyt/] acting on End(V), where the tensor product V = (8> ■ • • C?) V/ represents the 



'"Compare, however, eq. (|68] | with the operator k*^"^ defined in Lemma A. 2 of [8]. 
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space of states of a segment of the infinite spin chain reaching from site k to site /, and Vj is 
isomorphic to C^. Then we use that these endomorphisms have a reduction property similar 
to (|3TI) which allows us to extend their action to Wa by an inductive limit procedure. The 
endomorphisms bj^/j, C[yt/] and H[;;./] ^re constructed from weighted traces of the elements of 
certain monodromy matrices related to Uq{sl2). These monodromy matrices are obtained from 
products of L-matrices with different auxiliary spaces. 

The simplest case is directly related to the i?-matrix of the six- vertex model, 



/l 














1/ 



(44) 



where 



1-^2^2 



Y(C) 



1 



(45) 



and q = e"^. Let us fix an auxiliary space Va isomorphic to C^. Then Laj{Q = RajiQ is the 
standard L-matrix of the six- vertex model. The corresponding monodromy matrix is 



(46) 



It acts on Va ® V. We are interested in operators acting on EndCV). Such type of operators are 
naturally given by the adjoint action of operators acting on V. An example is the transfer matrix 
t[yt,/](0 defined by 

t[,,,](0(X) =tr„r,,[,,;](0-ixr,,[,,,](0 (47) 

for all X E EndCV). It will be needed in the definition of the operator ;] below. 

Further following [6] we introduce another type of monodromy matrices for which the aux- 
iliary space is replaced with the ^-oscillator algebra Osc generated by a,a*,q'^^ modulo the 
relations 



D * 
q a 



a a 



1 



-2D 



q^a 



aa 



aq 
1 



2Z)H 



We consider two irreducible modules W"^ of Osc, 

W+ = ff)C\k), w- 



k>0 



^C\k), 
k<-l 



defined by the action 

q^\k)=q'\k). 



a\k) 



J-k\ 



\k-\), a*\k) = {\-h,-i)\k+\) 



of the generators. The L-operators L {Q G Osc®^nd{V) are defined by 



L 
L 



(0 
(0 



•^-1/2^-1/4(1 _^^*^H 



C,ao - C, 



2J2D+2 



(5-0+)q' 



(48) 
(49) 

(50) 

(51a) 
(51b) 
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The corresponding monodromy matrices are 

where the index A refers to the auxiliary space Osc. We denote their (inverse) adjoint action by 

^tm^Q-'i^) = T^MiQ-'xT^^MiQ (53) 

for all X E End(V). Here the inverse on the right hand side is taken for both auxiliary and 
'quantum' space. The analogue of the transfer matrix /] in this case are two Q-operators 2^ 
(see [6]). Since we need only one of them here we leave out the superscript and defined 

Q[,,,](C,a)=tr+(^2aD,T+^^^^^(^)-i)_ (54) 

Here tr^ signifies that the trace is taken over W^. Similarly we will denote the trace over W 
bytr^. 

Now we are prepared to define the restriction of the operator H to End(V), 

%,/](C;a) = Q[,,,](C;a)t[,,z](0. (55) 

We show below that this definition (in the limit a 0) can be inductively extended to Wa- To 
avoid possible confusion let us note that in fact the operator H defined by the formula (l55l) is not 
the left hand side of Baxter's TQ-relation. In order that it were we would need to 'a-deform' 
the t-operator as well. 

In order to obtain bj^;. /] and /] and also another form of the operator /] we recall the 
fusion technique used in [6]. There the fused L-operators 

= iG%XKj^QRa,j{QG%^ (56) 

were defined, where 

The application of G^^^ transforms j{QRaj{Q into a matrix of lower triangular form on Va, 



The inverse is also of lower triangular form and is given by 



llHere we use a slightly different definition of Q-operator in comparison with Q+ in [6], see formula (2.10) 
there. The difference is an additional factor (1 — <7^'"^^'). 
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Correspondingly 

is of upper triangular form. It follows that similar statements hold for the monodromy matrices 

= (^t)"'^At[M^^)^«>M(Q^t • (61) 

^{A a} [k l\ ^'^^^ ^ lower triangular matrix in Va, (Q as an upper triangular ma- 

trix. As before we are interested in the adjoint action of the fused monodromy matrices on 
endomorphisms Z e End(V). Following [6] we define 

forallZ eEnd(V). 

Regarding [A: /] (^)~^ matrices acting on Va as in [6] we may write their entries as 

-^cx,,„(o k,jp): 

The entries of these matrices are elements of Osc ®End{V). We are now prepared to define 
^[k,i] and C[fc^;], 

C[M] (C, cx) = (1 - sing tr+ (?^«^^C+ (C))] , (64a) 

b[.,/](C,a) = ^''l^-'l sing[C-«+V/ltr-(^-2«(^-+i)B-j,^^j(Q)] . (64b) 

The symbol 'sing' means taking the singular part at ^ = ^j, 7 = 1, . . . ,n (cf. eq. (2.13) of [6]). 
These operators raise or lower the z-component of the total spin by one, 

[§M,C[;t,/](C,a)] = C[fe,/](C,a), [§[^,;],b[^,;](C,a)] = -h[k,i]i^,a) . (65) 

Their properties were extensively studied in [6, 8]. Here we shall only need the following. 

Proposition 1. Reduction properties [6]. 

b[k,i] (C, a) (%/-!] /;) = b[;t,z-i] (C, cc) // , 

CM(C,a)(?""^*%+i,/]) = ?(«-^)^*C[,+i,,](C,a)(Z[,+i,,]) , 
b[A,/] (C, «) = b[,+i,,] (C, a) {X^k+i,i]) ■ (66) 
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From this it follows that C[;;.,/](CiOc) be inductively extended to an operator c(^, a) : 
Wa W(x_i. Similarly b[^;](^,a) inductively extends to an operator b(^,a) : Wa Wa+i. 
These are the operators appearing in the definition (I38al) of Qi. 

Using the simple relation 

(GiJ-V'^'Gl^ = ,"""{l ('-^"'>3)^ (67) 

and the concrete form of L^^ ^| j{Q and L^^ '■^^^ obtain 

H[,,,](C;a) ^ {l-q-'^)t4{q^^-aXC+^^,^{Q) , (68) 

where the symbol ^ means equality up to the regular part when C, ^ ^j. Since the function 
(p(/j, a) is regular when fj ^ 0, the regular part of H[^ ;](^;a) does not contribute to the right 
hand side of (I38bl) . The formula (|68l) looks rather similar to the definition (I64al) of the operator 
c^fc.i]- The essential difference is due to the insertion of a\ under the trace. In contrast to the 
C[yt,/] -operator which increases the total spin, the operator tl[k,i] does not change the total spin. 



Properties of the operators Q.\ and Q.2 

Assuming for a moment that the limit on the right hand side of (I38al) exists we can conclude 
with (I66l) that 

{^i)[kj]{X[k,i-i]Ii) = {^i)[k,i-i] 

(^l)[A:,/](4%+l,/]) = h{^i)[k+\,i]iX[k+i,i]) ■ (69) 
Due to this property one can define D.i as the inductive limit of its restriction 

= lim lim(^2i)r^,i. (70) 

As we shall discuss later the same is also true for the operator ^2- 

But before we come to this point let us check whether the limits in the right hand side of 
(|38al) and (|38bl) are really well defined. 

Proposition 2. The limits in the right hand side ofeqs. ^38a^ and ^38b^ exist. 

Proof. The existence of the limit in (|38bl) follows from the formula (l68l) . because taking the 
trace there can results in at most a simple pole 1/(1 — q'^)- This pole will be canceled by the 
factor (1 — q^^'^) which stands in front of the trace in (l68l) . 

In order to prove the existence of the limit in (I38al) we use an alternative representation 

' ://|^|^(|)"<o(...,;a,XK„foa 



Q.\ = — lim 



(71) 
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wheretl 

X(^i,C2;a) = sing^^,^^[tv,(5,,,(^i/^2)T,(C2)-^T,(^i)-i)Q-(^2;a)Q+(Ci;a)] (72) 
with the 'boundary' matrix 



B{Q 



/o 








o\ 







-q-q^^ 










C + r^ 













0/ 



(73) 



and the same operators as defined in [6], 



(74a) 
(74b) 



The form (ItTI) of is similar to the form shown in the appendix of [6]. It can be obtained 
combining the ideas of [4] and [6] . 

The limit in (ItTI) exists, since the integrand is antisymmetric in , ^2 in the limit a ^ 0. This 
can be seen as follows. First of all (0(^25/^1 JOc) (^1/^2)" is symmetric in ^1, ^2 for a ^ (see 
eq. (fT9l)). Next tr^j, {Ba,b{^\/^i)'^b{^i)^^'^a{^i)^^) is independent of a and antisymmetric in 
^1, t,2, since 5(^1/^2) is antisymmetric in ^1, ^2 and since [5(^i),_/?(^2)] = 0- 

It remains to show that Q^(^2;oc)Q^(Ci;oc) is symmetric for a — 0. This product is mero- 
morphic in a by construction. We show by an explicit calcualtion in appendix |A]that it is regular 
at a = and symmetric in ^1, ^2 in this point. In fact, adopting the notation 



0[ ,...,CT„;e[ ,...,£„ 

for the matrix elements of the operators Q^(^;0) with respect to the canonical basis we obtain 

^\l.l<j<k<n{{^i-^'jK-i^j-^j)<^'k) . (76) 



Hence, 

Q-(C2;0)Q+(Ci;0) 

= Q+(C2;0)Q+(Ci;0) = Q+(Ci;0)Q+(C2;0)=Q-(Ci;0)Q+(C2;0), (77) 

where we used the commutativity [Q+(^i;a), Q+(^2;oc)] = (see [8]) in the second equation. 

□ 



*Here we take only the spin-0 sector. 
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Following the same lines one can show that the operator Qi is symmetric under the spin 
reversal transformation, 

ai = miS. (78) 

Moreover, Q.i is symmetric under reversal of the direction of the magnetic field 

^1=^4^-.' (79) 
since (O is an even function of the magnetic field h. An actual calculation of the right hand 
side of eq. (|38a|) or dTTI ) demands to apply I'Hopital's rule. As a result one gets two terms: 
one standing with (0(/Ji,/J2;0) which is even with respect to the transposition of fji and fj2 and 
another one with (d'{iji,ij2',0) which is odd with respect to /Ji ^ ^2- This is the same splitting 
as discussed in the paper [7]. Below in section [5] we will consider several examples in order to 
illustrate this point. 

Let us now come to the properties of the operator O.2. We shall consider 

Hy(^i'a\ ■ • • enl",) = res^.^^. H[i,„] {e,l\ . . . . (80) 
In the following we shall need an explicit formula which is also proved in appendix lAl 

Hi(ei^;...e4;J = (Q^\Mi;2,.,„)(e2^2,...e4;,), (81) 
where the action of the operator Mi;2,--,n is defined by 

IRl;2,--,n(^[2,)i]) = Rl,! ■ ■ ■Rn,\^[2,n]Rl,n ■ ■ -Rl,! (82) 

with the standard i?-matrix of the six-vertex model Ri j = Rij{b,i/b,j) and where the matrix 
elements of the operator Q^j are explicitly given by 

L*^OiJo'2,...,o;,;e^,...,e;,-^ei+---+e"'^i+---+^«%+---+e;,,<^2+---+<^« 



X 



Note that the limit a — > and the calculation of the residue at ^ = in equation (ISTI) may not 
be interchanged. 

The a = limit of the residues at ^ = for j >2 can be obtained from the formula (fSTT ) by 
applying the exchange relations 

Mu+iH[i,„](C;a)(X[i,„])=Hf;;;+i)(C;a)M,-,-+i(X[^^ (84) 
with H[!''t = Hri „i L . and the action 

= Ri,i+lX[l,n]Ri,i+l (85) 

for 1 < z, z + 1 <n. For example, 

H2(ei^, . ..e^l) = I^i^MF^^Xki/^ .^) {ei%e,% . ..e^^) . (86) 

A most important consequence of the explicit formula (l83l) is the reduction property 
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Proposition 3. 

Hi{he2l\...enll)=0, (87a) 

H, (/i e2l\ . . . enl) = h H, {e2% • • • ) , 2 < j < n , (87b) 

H; {e,l\... en-i%";\ In) = H, {e,%\ . . . en-i%\-\ ) In, 1 < 7 < n - 1 , (87c) 

lln{eil\...en-il:;_\ln)=0. (87d) 



Proof. The first formula (I87al) is rather trivial because from the formula ([83]) it follows that 

ETn^l ^2, •••,On;e2, •••,£« _ D 

CJ=±1 

The second formula (|87b|) is less trivial. Let us outline the proof for j — 2. First we use (|86l) in 
order to obtain 

H2 (/l .2^ . . . = ^2) (q| M2;3,...,„) (^4^31 . . . en%,) (88) 

and substitute eq. (f83l) . The latter should be separated into two parts in such a way that only one 
of them is touched by two i?-matrices in the right hand side of (f88l) . This part looks like 

where = i^^n and where the indices £3, . . . ,e', are considered to be fixed. The following 

identity can be verified directly, for example, on a computer 

If we substitute the right hand side back into (l88l) and collect all pieces we come to the statement 
that we wanted to prove, namely, 

H2(/i e2|...e,4) =hll2{e2l\...enl). 

The other cases when j > 2 can be treated in a similar way. The formulae (I87cl) . (|87dl) are 
simple consequences of the inversion of L-operators in the definition (l55l) . □ 

Using proposition [3] one immediately comes to the reduction relation for ^^2 because the 
restriction of (I38bl) to the interval [l,n] is 

(^2)[i,„] = -l9(^;;0)H,. (89) 
7=1 
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Proposition 4. Reduction identity for Q.2. 

(^2)[l,„] In) = (f^2)[L„-l] (X[l,„-1]) In, 

{a2)[Un] {hX[2,n]) =h {^2)[2,n] (%,„]) " (90) 

Due to dPOl) we may define ^2 for the infinite chain through an inductive limit as in eq. dTOl) . 

Another immediate consequence of the formula (|83l) is the spin reversal anti-symmetry. First 
of all 

L^-oiJ-o^,... -o;,;-e^,...,-£;- L^f^iJa^,...,a;,;e^,...,e;, " '^^^^ 
Then, since the operator Mi;2,...,n is symmetric with respect to the spin reversal transformation, 

IRl;2,...,« = J[2,n]IRl;2,...,«J[2,n] , (92) 

the operator Hi defined by (fSTI) is spin reversal anti-symmetric 

Hi = -JHiJ. (93) 
The same is true for the other residues with j > 2. Hence, one concludes that 

Q.2 = -J^2J- (94) 

Moreover, due to the fact that the function cp given by eq. ([T?] ) is an odd function of the magnetic 
field we have 

^2 = -^2|,_,. (95) 

The splitting of the whole operator Q in equation (|37] ) into two terms Q.\ and Q.2 seems rather 
natural because the two terms are even and odd with respect to the reversal of the spin and the 
magnetic field, respectively. 

5 Examples 

Li this section we present explicit formulae for the density matrices for n = 1 , 2 and for some 
particular matrix elements and correlation functions for n = 3. Since the definition of the op- 
erators b, c and H involves the multiplication of In two-by-two matrices and subsequently the 
calculation of the traces over or , it is already cumbersome to work out by hand the case 
n — 2. We preferred to use a little computer algebra programme for this task. 

The case n = 1 

This case is rather simple because f^i = and Q = Q.2. Since £i} = Q^ = Q one should expand 
the exponent in eq. (l36l) only up to the first order with respect to O.. A direct calculation shows 
that the operator Hi acts on the basis elements as follows. 



Hi(ei±)=±^/i, Hi(ei^)=0. 



(96) 
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Then from (|89| ) one obtains by multiplying the above result by — cp(?ii;0). It is left to 
substitute it into the formula (l36l) and take the trace ^tri . Finally one obtains the inhomogeneous 
density matrix 

D,(X,|r,A;0) = i/,-!?i^a^ (97) 

In particular, setting Xi =0 one obtains (see ([5]), (l8al) ) for (twice) the magnetization 

(o^)r,/, = tri(Z)i(r,/i)o^i) =-(p(0;0). (98) 
This result is in full agreement with equation (74) of [22]. 

The case n = 2 

This case is already less trivial. First let us calculate Using I'Hopital's rule and the fact that 
the functions (o(/ji,/J2;0) and (jo'(/ji,/72;0) (recall the definition (l20l) of oo'!) are even and odd, 
respectively, with respect to the transposition of fji and /J2 (see eq. (fT9l) ) one obtains 



where 



= -(0(Xi , X2; 0) a+ - co'(Xi , X2; 0) , (99) 



a+ = Inn (bi(a- l)c2(a) (^i/^2)" + b2(a- l)ci(a) feAi)'^) , 

= lim a(bi (a - l)c2(a) (^1/^2)" - b2(a - l)ci (a) (^2/^1)'') , (100) 



and 



b;(a) =res^^^.|^b(^;a)^^ , c^a) = res^^^. |^c(^;a)^^ . (101) 
The result of applying the operators Q.f to the basis of the 5" = sector is 

a+(.il.2g) =-f cth(Ti)A/2, ^i{e,_le2-D = \'-^^^^hh, 
£o 1 chfn) 

^^r(^le^2S) =-^Cth(Xi-X2)/l/2, ^li^l-le2-D=^ ^^^^^^\^^ hh. (102) 



It is clear that 



{aff = a+a- = a-a^ ^0 (i03) 



which implies 

Also the symmetry with respect to spin reversal is obvious in the above explicit formulae (11021) . 



Q.l = 0. (104) 
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Let us proceed with the anti- symmetric part. To obtain Hj for j = 1 , 2 one can either take 
the corresponding residues in the formula (l68l) or one can use the formulae (ISTI) for j = 1 and 
(|86l) for j = 2. The result is 

Hi(.i|.2^) = |(/r(^i,^2).ii.2i+/r"(^i, ^2)^1:1^2:1+ 

+/f(^i,^2)^^iIe2:|+/2-'"(^i,^2)ei:|e2l-a^i(^i,^2)(ei;^^2i-eil^^2;)), 
Hi(.r|.24) = ^(^-VK^i, ^2)^11.21+^/3-^(^1, ^2)^1:1^2:1+ 

+ £^"^gj(^i,^2)ei|e2:|-e^^2 (^I'^2)ei:|e2| + ,?3(^i,^2)(ei^e2i-eile2+)) (105) 

and 

H2Ke2g)=|(/f(^i,^2)ei|/2 + /+(^i,^2)ei:|/2), 

H2(er|e24)=-|(^-V3+(^i,^2)eiI/2 + ^/3-(^i,^2)ei:|/2), (106) 



where 



/i+(^l,^2) 
/2+(^l,^2) 
/3+(^l,^2) 



l_^2/^2 ' 

(^-^-1)^ + (1-^1A1)^ 

(l-^2/^2)^^^^/^^_^-l^^/^^)(^^^/^^_^-l^^/^2) 

1 



and /r(^i,^2):=/-+(^2,^i). (107) 



^1(^1,^: 



2 



(^iA2+w^i)(^-^-^: 



(^^1/^2 - ^-iW^i) (^W^i - ^-1^1 A2) ' 

g . _ (^-^-^)^+g^^(^iA2-w^i)^ 

^2l<;i.<;2j. (^^/^2-W^l)(^^lA2-^-lW^l)(^W^l-^-l^lA2) 



2 —2 

^'^^^'^'^ ^= (^^lA2-^-lW^l)(^W^l-^-l^lA2) ■ ^^^^^ 

The anti-symmetry of the operators Hi and H2 with respect to the spin reversal transformation 
is evident in the above formulae. 
Also one can directly verify that 

Hf = H^ = HiH2 + H2Hi =0 (109) 

and 

Hj^li+aiHj = 0, j = l,2. (110) 
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This means that the operator Q2 which is 

^^2 = -cp(^i;0)Hi-(p(X2;0)H2 

satisfies 



From this follows that 



^^2 = 



(111) 

(112) 
(113) 



and the expansion of the exponent in the formula (|36l ) extends only up to the first order in 
powers of Q.. 

Therefore in order to compute the elements of the density matrix we need to calculate the 
traces 

^ ^ - ' (114) 



Di l\% iliM\T,h;0) = - tn tr2 [(id (eie^/ e2^e^)] . 



For this purpose we have to use the formulae dlH), (11021) and (IlllL (11051) . (11061) . The result 
decomposes as follows, 



D2{h,X2\T,h;0)=Df''^{hM)+Df\hM) 



(115) 



where and ^2'^'^ are 4 x 4 matrices, 



Dr"(XiA2) = J/®/+^ 



ch(Xi-X2) .^^ , ch(Tl) ^j^^ . 

(o(Ai, A2;0) H — —00 (Ai, A2;0) 



sh(ri) 



rish(?ii - ?i2) 



(a+(g)a +a 00+) (116) 



and 

^r(^iA2) 



,odd,....^_ ^(^i'Q)aZQ/ ^(^2;0) ,^^, 
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sh(ri)((p(Xi;0)-(p(?i2;0)) 



(a+(g)0 -a (g)0+). (117) 



4sh(?ii -?i2) 

The homogeneous limit Xi,X2-^0 can be readily taken. We obtain the density matrix for n = 2. 



D2{T,h) = - 



/(g)/-(p(o^(g)/+/(g)0^) -sh(ri)(p;,(o+®o -o (g)0+) 

00 ch(ri)(0^ 



+ ( cth(ri)a)+ ) a^(g)o 



sh(ri) 



+ ■ 



where we introduced the shorthand notation 



(p = (p(0;0), cp.^ = axcp(?i;0) , (o = (o(0,0;0) , d^ = dx,(£l {XiM\^) 

A=0 



(118) 



• (119) 
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The density matrix (II 181) can now be used to obtain any two-site correlation function, e.g., 



{o\ol)T,h = tri2{D2{T,h)a\al)=cth{r[)(o+^, (120a) 
((T{a^)r,. = tri2(D2(r,/.)(T{al) =-2^-^^^. (120b) 

The case n = 3 

The explicit forms of Q.j or Hj are already quite involved for n = 3. We shall not present the 
exhausting list of matrix elements, but rather restrict ourselves to some examples of physical 
interest. 

We introduce shorthand notations 

f?''^ = /fn^2,^3)/rfe,^i) , ^i"'' = /r (^i,^2)/rfe,^i) , 

t^"^^'^^ =/fn^l,^2)/f^fe,^3)/r(^3,^l). (121) 

Using these symbols, the longitudinal correlation is represented rather compactly. In the inho- 
mogeneous case we find 

= th(ri) + op) a)(Xi M;0)+ th(ri) + ©(Xz, h;0) 
+ (2 cth(2ri) + c)^") + cth(ri) (0+^ + O3 +)) a)(?i3 Ai ; 0) 
-t+++(a)'(^iA2;0)+a)'(^2A3;0)) 



4sh2(ri) 



( 



- - (4ch2(ri)t+++ - (t++- + {-+- + 



+ t+-+ + r++ + r-+)(o'(;i3,?ti;0). (122) 

Taking the homogeneous limit we arrive at 



lid)' sh^fn) 
{o\ol)T,h = 2cth(2ri)a)+ -th(ri)a),,- -th(ri)a).^ + ^ - (123) 



By X and y we denote the derivatives with respect to first and second argument taken at zero. 
The same limit for the transverse correlation reads as follows. 



1 , ch(2ri)th(ri) ^^ ch(2ri)th(ri) ^^ ch(2ri) ^^, , sh2(ri) ^^, 

:^ -4 ^xy 7; a>xx ^r— % H — ?r— a>xxv ■ (124) 



sh(2Ti)"^^ 4 "^^^ 8 2ri 8ri 
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The rational limit in the last two equations is not easy. Using the high-temperature expansion 
we checked to 0( l/T) that it coincides with our previous result [1] for the XXX chain. 

As a last example we show the emptiness formation probability in the inhomogeneous case, 

D+++(XiA2A3|r,/i;0) = ^ + ^(-(p(X3;0)+Ci(^i,^2,^3)«(^i,^2;0) 

+ C2(^l,^2,^3)w'(^lA2;0)+C3(^l,^2,^3)«(^lA2;0)cp(^3;0) 

+ C4(^i,^2,^3)(o'(?^i A2;0)cp(?i3;0) + cyclic permutations ) . (125) 
Here the coefficients are given as follows, 

Ci(^i,^2,^3) = (2cth(2ri)(0++ + c)r)+cth(ri)(0+-+fr+)), 

C2(^i , ^2, ^3) = ^ (-2ch(2ri) t+++ - (t++- + + 1+-) + 1-+ + r +- + r++) , 

C3(^i,^2,^3) = cth(ri) (2(0++ +Dr") - (j^r +f r^)) . 

C4(^i , ^2, ^3) = ^ (-4ch2(n) t+++ + 1++- + 1+-+ + r++ - (t-+ + r+- + 1+-)) . (i26) 

The homogeneous limit is left as an exercise to the reader. 

6 Conclusions 

In an attempt to generalize the recent results [6, 8] on the factorization of the ground state 
correlation functions of the XXZ chain to include finite temperatures and a finite longitudinal 
magnetic field we have constructed a conjectural exponential formula (|37l) . (|38l) for the density 
matrix. The main steps in our work were the construction of the operator H, eq. (|55l) . which 
takes care of the modification of the algebraic part of the exponential formula in the presence 
of a magnetic field, and of the functions cp and (O, eqs. (fTSi) . (fTTI) . which allowed us to give a 
description of the physical part in close analogy to [6, 8]. In the limit T,h ^ our conjecture 
reduces to the result of [6, 8], even for finite a. It also trivializes in the expected way as T ^ oo. 
We tested our conjecture against the multiple integral formula ^ by direct comparison for 
n = 2 (see appendix |B]) and by comparison of the high temperature expansion data for n = 3 
and n = 4. Judging from our experience with the isotropic case [17] we find it likely that very 
similar formulae also hold in the finite length case and that the only modifications necessary 
to cover this case are a restriction of Q. to the finite length L of the chain and a change of the 
auxiliary function from (flOl) to (fT4)) . 
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Here we outline the proof of eqs. (1761) and (l83l) . Our starting point is eq. (|74al) . Since we work 
in the sector = 0, we have to set S = 0. Then 



where (L^(Q)^ are the matrix elements of the L-operator (I51al) , 



e,e' 



and 

The main observation is that for the computation of the limit a ^ of eq. (lA.ll) it is enough to 
substitute there and (L^y^ by and {L^y^ with(lll 

and 

where we set 

aj" = aA , a^=a\. (A. 6) 

In this notation the algebra (l48l) looks very simple 

The reason is as follows. Let us first formally substitute and (^^)^^ for and (^^)^^ into 
the right hand side of eq. (|A.1I) . 



_ 1 

X 



^^Strictly speaking the operators and (L^) ^ are not inverse to each other any longer. 



26 



We do not write the 5's like in the right hand side of (|76|) which reflect the fact that we are in 
the spin-0 sector. Let us just imply that they are there. 

Let us formally ignore all inside the trace here. Then the total degree of q^'^ is zero 
because £y e'j = Y,j and only q^^^^ i^fi^ which produces a term 1/(1— q^^) after taking the 
trace over the oscillator space A. Since the differences between and and between {L'^)~^ 
and (L^)^^ contain only positive powers of q^-^, the insertion of such terms does not change 
that most singular term 1/(1 — q^^) when a — 0. Therefore we can ignore those differences 
when calculating the limit a ^ 0. 

One more observation is about the contribution coming from the terms containing a^. Sup- 
pose we had just 

tr+(4'...al>2'^-) 

with £i H \-e2n =0. Then, using the algebra (l48l) we would conclude that again the most 

singular term would be 1/(1— q^"') as a result of taking the trace. It means that if one succeeds 
in collecting all then one can replace them by 1 without any change in the most singular 
term. The first conclusion obtained from the above is that the limit a ^ of the expression 
(IA.8I) gives us the limit a — of eq. (lA.ll ). Second, in order to calculate it we have to collect all 

inside the trace (IA.8I) using the algebra (IA.7I ) in one place, say in the place of the symbol x in 
(IA.8I) . If we do this and afterwards ignore the product of all following the above arguments, 
then we can easily take the limit a ^ and come to the formula (1761) . Similar arguments may 
be applied when treating the a ^ limit of the formula (|74b|) . 

Now we outline the derivation of the formula (|83l) . When calculating the residue at which 
is implied in eq. (|8T1) one obtains 

(L+(CA,))^;...(L+((;A,))^S(i+(C/5,))°r°"') , (A.9) 

where summation over a is implied. The pole at ^ = originates from the L-operators with 
argument ^/^i. We use the cyclicity of the trace and directly verify that 

res^.^,(L+(C/^0)^^"^^(4(C/^i)-^):^ 

= ^[^-aj^^(^ ^^'^ ) 

Implying that we need to calculate the limit a ^ in the end we may set 

^2«Oa = Q _^a^£^^-i(ei-ai)^2aD^ ^^^^^ 

ei,Oi 



^laD^ (A. 10) 



1 



(1.4) 



on the right hand side of (lA.lOl) . Thus, we come to the conclusion that the right hand side of 
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(IA.9I) is equal to 
lim 



(i+(5,/5,))'"...(L+K,/fc))';e,a-i<='-<"V"''')l . (A.12) 



Finally we can apply to eq. (|A.12I) the same trick as described above in order to get the for- 
mula (1831). 



Appendix B: Factorization of the double integral 

In this appendix we show that our conjectured formula for the density matrix for n = 2, eqs. 
(|115l )- (|117| ). coincides with the double integral, eq. Q for n = 2. The density matrix for n = 2 
has six non- vanishing elements. In this appendix we will denote it by D rather than by D2 and 
suppress the temperature, magnetic field and a dependence of the matrix elements for short. 
Using the Yang-Baxter algebra and reduction we find four independent relations between the 
six non-vanishing matrix elements of D, 

D+ZiliM)=D+ili)-D+X{liM). DzXiliM)=DXil2)-D++{li,l2), 
DZZihM) = D++{liM)-D+ih) -DXil2) + 1 , 

D;l(.a.)-ol;(.a.)^«;^i^iM. (B.I) 

sh(Ai— A2) 

Inserting these relations into ?i2|r, h;0) we obtain 

Z)(Xi,X2|r,/z;0) = ^/®/+^(2D+(;ii)-l)o~®/-f^(2D+(X2)-l)/®o^ 
+ ^ {4DXX ihM)~2DXili)-2DX{l2) + l)o'0o' 

+ 2,^x^1 X,) (^1= ^^1) - ^^2)) (a+ ® a- - a- ® a+) , (B.2) 

and we are left with the problem of expressing the one-point function and the two-point 

functions DXX (^1 A2) and (Xi , X2) + (ki , ?i2) in terms of cp and co. 

Comparing dH) for n = I with the definition (fTSi) of our function cp we find the relation 

2DX{X) = l-(?{kO) (B.3) 

for the one-point function. 
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In order to simplify our task for the two-point functions we introduce the quantum group 
invariant combination [4] 

= e^i -^2^1; {Xi ,X2)+ e^^-^' D-t (Xi A2) - e-^D+I (Xi A2) - e^Z)I+ (Xi A2) • (B.4) 
Using again (|B.1I) we obtain the relation 

= DqihM)+cH^){DXi'kl)+D+{l2)-2D++{hM)). (B.5) 

Hence, in order to determine the density matrix for n = 2, it suffices to calculate Dq{X\ , X2) and 
D^^(?ii,X2) from the double integrals. 

Let us start with the simpler case Dq{X\,X2)- Inserting ^ into the definition (IB .41) we find 

Dq{XiM)= [ / ^ det[-G(a);,;i/^;0)]KtOi,tO2), (B.6) 

7e 27ii(l + o((0i)) Je 2711(1 + 0(0)2)) l ^ ^ /j ^ 

where 

e^'+^^sh(Xi-X2+Tl)sh(Xi-X2-Tl) 

K»i,»2) = - eco,+co.,h(;.i-;.2)sh(a)i -(02-11) • ^^-'^ 



Using the simple relation 
we can rewrite (|B.6I) as 



^ + . L ^ = 1 (B.8) 



l + a(a)) l + a(a)) 



/),(^i,^2) = ^ 27ti(lTa(a))) (Kto,Xi)G(a),X2;0) -r(co,X2)G(a),Xi;0)) 



'e 27ri(l + a(a))) 
1 /" d(Oi f da)2 



■det[-G((O^Afc;0)] (r((0i,(02) -r(0)2,G)i)) . (B.9) 



2Je 2ra(l + a((0i)) 27ri(l + 0(002)) 

The first term on the right hand side is already a single integral. For the second term we observe 
that 

1 ch(Ti)(e-2"i -e'2co2) sh{h-l2+y])sHh -l2-y]) 

2^ ^ (e-2A.i _e-2'^2)sh((Oi-(02+ri)sh((Oi-(02-ri) 

The (D-dependent terms in the denominator are proportional to the kernel in the integral equation 
(fTT)) for a = 0, and the numerator is a sum of a function of (Oi and a function of a)2- Hence, 
the double integral can be reduced to single integrals by means of the integral equation (fTT)) . 
Collecting the resulting terms and inserting the definition (fT6l) of our function \\f we arrive at 



sh(?ii - X2 + ri) sh(?ii - ?i2 - ri) ^ 
2sh(ri) 



Dq{hM)= A2; 0) . (B. 1 1) 
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Let us proceed with the calculation of D~l'^{'ki,'k2) which according to ^ is equal to 

dcoi e-«^ r d(02 e-«^ 



e 27ri(l + a((Oi)) 27ri(l + 0(002)) 

det[-GKa>;a)] tir'>Vj)'"'^"''"t <B-'2) 
sh(A2 - Ai) sh((Oi -(O2 -Tl) 

^ V ' 

=i-(C0i,C02) 

Because of the antisymmetry of the determinant we may replace ^((Oi, (O2) with 

^(5(0)1,(02) -^(a)2,(0i)) 

ch(ri)(sh(2a)2-?ti -'^2-'^) -sh(2(0i -Xi -?i2-Tl)) +ch(?ii -?i2) sh(2(a)i -002)) 
4sh(?ii -?i2)sh((0i -(O2 + T1) sh((Oi -0)2 -T]) 

Then 

^'++ 

where 



D+X A2) = ^1 A2) + lim 72 A2; a) , (B. 13) 



7l(?ll,?l2) 



72(?^i,?^2;oc) 



dCO; 



2 ^ 

nye27ri(l + a(a);)) 

r 2 . do), e-«^ 
IJie 27ri(l + a((0,)) 



det [G{iOj, Xk; 0)] ch(ri) sh(2(02 -Xi-X2-r[) 
2 sh((Oi - (02 + Tl) sh((Oi - (02 - Tj) sh(?ii - X2) 

det[G(a)j,?iyt;cx)] cth(?ii -X2) sh(2(a)i -0)2)) 
4sh((0i -a)2 + Ti)sh((0i -(O2-T1) 



Here 7i {Xi , ?i2) is of a form which allows us to carry out one integration by means of ([TT]) (for 
a = 0). The result is 

^1(^1,^2) = \ - ^(cp(Xi;0) + (p(X2;0)) - ^cth(ri)i|/(^iA2;0) 

1 , '^ N v-i . /■ da)G(a)A/>i;0)cth((o-?ip2-T|) , 
+ -cth(Xi-X2) £ sign(/') / " ^^-l"^) 

2 ^1^2 7e 2711(1 + o((o)) 

For the calculation of 72 (Xi , X2 ; oc) we express the hyperbolic functions in the integrand in terms 
of the kernel (fT2l) occurring in the integral equation (fTT)) for G, 

sh(2((Di -(D2)) ^ i^((oi-(D2;«) -i^((02-(0i;a) 

sh((Oi -(02+Ti)sh((0i -(O2-T1) 2sh(ari) 

Then the integral over (O2 can be performed by means of the integral equation ([TT]) for finite a, 
and we obtain 

4 2sh(ari) 

-ie--''cth(^,-X,) £ sig„(P) / dooG(co.V.;a)c.h(oo-V.-,) 
2 Je 2711(1 + a((o)) 
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From the definition (fT6l) of \|/ and from the integral equation (fTTI) for G we infer the symmetry 
property \|/(?i2,^i;oc) = V(^iA2;— oc) which can be used to carry out the limit a — for J2. 
Using it and inserting the a ^ limit of (|B.16I) and (IB. 141) into (IB.13I) we arrive at 

= \ - i(cp(^i;0) + cp(^2;0)) - ^V(^a2;0) - ''^^^^M^\X,M;Q) , (B.17) 

where the prime denotes the derivative with respect to a. Inserting now (IB.5I) . (IB.l II) and (IB.17I) 
into (IB .21) and taking into account the definitions (flTl) and (|20l) of (0 and (£>' the reader will readily 
reproduce the density matrix (|115I) - (I117I) for n = 2. 



Appendix C: The high temperature expansions 

We comment on the application of high temperature expansions (HTE) to the multiple integral 
formula, which provide important data for the construction of the conjectures in this report. This 
may also be a basis for the numerical evaluation of correlations as demonstrated in [41,42]. 

As is usual, we assume an expansion of quantities in regular powers of j. We then typically 
face the problem of solving a linear integral equation for a unknown function f{X), 

m=fQ{l)+v [ ^i^(X-(o;a)/(a)):=/o(X)+viS:*/(X), (C.l) 
Je 2%i 

where v stands for some constant. The driving term fo{X) is a known function which has at 
most simple poles at X = fjj and a pole of certain order at X = inside C. Eq. (|C.1I) can be 
solved in an iterative manner, 

m = fo{X) +vK*MX) Wk* {K*fo){X) + ■■■. (C.2) 

The crucial observation is that K* foCk) has poles at X = ±ri,/j/±ri and that these poles are 
outside of contour G. Thus, only the first two terms in (|C.2I) do not vanish and f(k) = fo{X) + 
vK*fo(k) solves eq. (IC.ll) . 

This mechanism makes it possible to evaluate each order in the HTE in an analytic and exact 
manner. Of course, the evaluation of residues becomes more and more involved with increasing 
order of j. Computer programs like Mathematica, however, can efficiently cope with such a 
task and we obtain sufficiently many data for our purpose. 

Here we present some examples which one can compute by hand. We consider the nonlinear 
integral equation (fTOl) under the assumption 

a(A) = 1 + -^r- + -^2- + ■■■■ 
Then comparing 0{y) terms, one obtains the equation, 

27sh2(Ti) 



aW(X)=ao{X)-l- [ ^if(X-a);0)a(i)(a)), ao{X) 
2 Je 2111 



sh(X) sh(;i + ri) ■ 
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We apply the above strategy and find the first thermal correction to a{X} as 



sh(X) sh(?i-ri)sh(X + ri) 



Similarly the first correction to a(A,) is found to be a^^\X) = —a^^\X). 

Equation ^ can be solved similarly. Let G(?i, /j; a) = G^^) (k, /j; a) + G^^^ /j; a) /r + ■ • 
Then the following explicit forms are obtained. 

G(«) {X, ,;a)=- coth(X e«^ ^-M^'^^-^) + cotHl-, + ^) 



+ 



4 2sh(^)sh(/j-ri)sh(/j + ri) 

7shrii^(;i;a)G(°)(0,A/;a) 



All elements of the density matrix can now be evaluated up to 0(r ^). A simple example 
is the emptiness formation probability for n = 2, 

n++n ^iTi,m 1 a(i)(Xi) + a(i)(X2) 7shriG(0)(0A2;0) sh^i sh(X2 -^i +Tl) 



4 8r 4T sh(;i2-?ti)sh(;ii-ri) 

_ /shriGW(OAi;0) shX2sh(Xi-X2 + ri) 2^ 
AT sh(Xi-X2)sh(?i2-ri) ^ 

The other basic functions are also readily evaluated. 

(P(a/; a) = - ^ + ^ (( 1 - e-'^^) coth(A/ - ti) + ( 1 - e«^) coth(A/ + ^))+0{T-'), 
1.. ^ (a(i)(^2)-a(i)(A/i))GW(A/2,w;0c) 



2 vr. r., . 2r 
- ^g(o) (0, ,.1 ; a) g(«) (^2, 0; a) + 0(r-2) . 

One can then check the validity of our conjecture by comparing the multiple integral formula 
for the density matrix and the exponential formula after substitution of the basic functions by 
their HTE data up to 0{T^^). The higher order terms can, in principle, be checked in the same 
manner. 

Before closing the paper, we sketch briefly how we used the HTE data to arrive at our con- 
jecture. Each density matrix element consists of two parts; D^^^^, the even part with respect to 
the magnetic field, and D°'^'^, the odd part. The factorization for n = 2 can be done fully in an an- 
alytic manner, as demonstrated in appendix |Bl This result and the previous results of the XXX 
case motivate us to assume that the even part shares the same algebraic part with the ground 
state case. Then it is not difficult to identify two basic functions, (|)(?ii,?i2)7<l^(^i,^2)- We can ac- 
tually represent them by the single function \|/(?ii,?i2;cx) such that (|)(?ii,?i2) = ^\(/(?ii,?i2;0). 
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(|)(A,i , X2) = — ^ ) ^2; cc) |a=o- This is one of the advantages in using the disorder param- 

eter a. 

We then consider the odd part of n = 3. The most interesting sector is Dg'^'^oV-oa^ (^1 ■, ^2, ^3) 
with — Y^i^i = 1 to which nine elements belong. With use of the Yang-Baxter relation and 
the intrinsic symmetry of the density matrix, one can represent all the element by only one 
element. We choose Dg'^'^l^^ for this, with permutations of the arguments (A,i,A,2, A,3). The 
resulting 6 objects are found to satisfy linear algebraic relations, and the consideration of the 
kernel space implies the representation 

nodd++-/e p p \ (^1,^2,^3) , 52(^1,^2,^3) 

^3 -++(qi,q2,q3) = ? H — ^ , 

S2 ';i';3 

where ^, = e^', and s\,S2 denote certain symmetric functions of 

We then assume that ^'1,52 are given by sums of products of rational functions of Xi and the 
basic functions (p, (|) and (|), e.g., 

51(^1,^2,^3) -Vo(^i,^2|^3)cpfe;0) + yi(^i,^2|^3)cpfe;OH(^i,^2) 

+ y\ (^1 , ^2 |^3)(p(^3; 0)^(^1 , ^2) + cyclic permutations. 

y/(^i,^2|^3) is symmetric in ^1,^2 0' = 0, 1) while yi(^i,^2|^3) is anti-symmetric. 

Furthermore, we restrict the possible forms of these coefficients according to our previous 
experience such that 

T/ e IP ^ P0(^l,^2|^3) „ .j: . \. x ;?i(^i,^2|^3) 

n)(qi,q2|q3) = 7^2 — pi\(Pi — E2V' n(qi,q2|q3)- 



Vl(^l,^2|^3) = 



(^f-^i)(^i-^i)' .v....^,.-^. (^2_^2)(^2_^2^ 

PI (^1,^21^3) 



(^?-^2)(^2_^2)(^2_^2)- 

Now pj and p\ are polynomials in ^( and pj (pi ) is symmetric (anti-symmetric) in ^1,^2- 
We then assume polynomials of certain orders with desired symmetry for them and fix the 
unknown coefficients so as to match the HTE data. All these parameters are fortunately fixed by 
the data up to 0{T^^). At the final stage, several hundreds of terms are canceled just by fixing 
one parameter, which looks rather convincing. We then check that the choice of parameters 
actually recovers the 0{T^^) terms in the HTE. After this procedure, we arrive at expressions 
for the density matrix elements now written in terms of rational functions and basic functions. 
We then try to fit them into the exponential formula. This requires the new operator H in the 
main body. Once it is identified, it is easy to write down the conjecture for n = 4. Then again we 
test the validity against the HTE data for the multiple integral formulae of the density matrix. 
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